Helium at elevated pressures: Quantum liquid with non-static shear rigidity 
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The properties of liquid helium have always been a fascinating subject to scientists. The phonon 
theory of liquids taking into account liquid non-static shear rigidity is employed here for studying 
internal energy and heat capacity of compressed liquid 4 He. We demonstrate good agreement of 
calculated and experimental heat capacity of liquid helium at elevated pressures and supercritical 
temperatures. Unexpectedly helium remains a quantum liquid at elevated pressures for a wide range 
of temperature supporting both longitudinal and transverse-like phonon excitations. We have found 
that in the very wide pressure range 5 MPa-500 MPa liquid helium near melting temperature is 
both solid-like and quantum. 
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INTRODUCTION 

This Letter is concerned with the thermodynamic func- 
tions of liquid 4 He and their possible relationships to 
those of its crystalline counterpart. Compared with crys- 
talline and gas phases the statistical description of the 
structure and thermal properties of liquids remains rel- 
atively incomplete. At low pressures, matter commonly 
exists either as a dense solid or as a dilute vapour. For 
each of these states there is a model which is a plausible 
approximation to reality and which constitutes a basis 
for detailed theoretical extension. These are the ideal 
periodic arrangements for crystalline structures and the 
ideal gas proposition, respectively: in the former, em- 
phasis is placed on structural order modified slightly by 
zero-point or thermal motion of the atoms while (save 
perhaps for the ultra-cold gases) the latter model de- 
scribes the thermal motion of the atoms on the basis of 
random atomic displacements and associated momenta. 
But as is well known there is a third state of dense matter, 
the ubiquitous liquid state drfl]. This state occurs over 
a temperature range that separates the regions occupied 
by the solid and vapour states. However, the problem of 
formulating a rigorous mathematical description of the 
molecular motions in liquids has always been regarded 
as much more difficult than that of the kinetic theory 
of gases or collective displacements of crystalline solids 
[7H16|. Approximation, if not judicious, can lead to a 
description of either high-density gases or of disordered 
high-temperature solids. Indeed, at one time consider- 
able effort was devoted to the representation of liquids in 
these terms, but it is now known that liquids do not have 
a simple interpolated status between gases and solids, 
although similarities to the properties of both adjacent 
phases can certainly be observed. 

The ability of liquids and solids to form free bounding 
surfaces obviously distinguish them from gases. The co- 



those of gases. And the viscosities of gases and liquids 
are some thirteen orders of magnitude lower than those of 
solids, and this we may easily understand in terms of the 
molecular processes of momentum exchange. In terms of 
vibrational states liquids differ from solids because they 
cannot support static shear stress. However, as will be 
seen below liquids support shear stress at high frequency. 
Flow in a solid arises primarily from rupturing of bonds 
and the propagation of dislocations and imperfections. In 
a liquid flow is characterized by both configurational and 
kinetic processes, whilst in a gas the flow is understood 
purely in terms of kinetic transport. In this very lim- 
ited sense liquids may have a minor partial interpolated 
status between gases and solids. 

The concept of elasticity and viscosity in liquids merit 
clarifications. Which property dominates, and what val- 
ues of the associated parameters are assumed, depends 
on the stress and duration of application of that stress. 
If we apply a stress over a very wide spectrum of time, or 
of frequency, we are able to observe liquid-like properties 



in solids and solid- like properties in liquids. Frenkel 17 1 
introduced a relaxation time r as the average time be- 
tween two consecutive local structural rearrangements in 
a liquid. If r is small compared with an observation time 
it will yield to the process of liquid flow. In this macro- 
scopic hydrodynamic picture, we now have a rather good 
understanding of most of the fundamental processes op- 
erating on such time-scales in liquids, including quantum 
liquids such as liquid helium. 

The investigation of the properties of helium has been 
one of the most prolific endeavors since its discovery (in 
1868) flil - Eoj ]. For most systems the solid phase is the 
state of lowest energy at one atmosphere. To date liquid 
helium is the sole exception in this respect: below 1.70° 
K liquid helium has a lower free energy than that of solid 
helium 
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efficients of self-diffusion of liquids ( 10 



) and 



22j • Liquid helium is a quantum liquid at low 
temperatures and quantum effects play a crucial role. He- 
lium is also of considerable practical significance which 



solids ( 10" 



) are orders of magnitude below is related to the rapidly growing industry surrounding 



2 



the various applications of superconductivity these often 
relying on liquid helium as a coolant or refrigerant. Con- 
siderable interest is still focused on the highly unusual 
properties of helium especially at low temperatures and 
the study of the properties of liquid 4 He continues to be 
an active area of condensed matter research 23142 7j . 

Theory has been long drawn to study the condensed 
isotopes of helium and their mixtures because these liq- 
uids are model many-body systems but with fundamental 
quantum-statistical differences; they are fertile proving 
grounds for various quantum many-body formalisms. As 
Landau emphasized [28[, these systems are amenable to 
theoretical attack because, when studied at relatively low 
temperatures, they are only weakly excited from their 
ground states. A description in terms of weakly interact- 
ing elementary excitations is then appropriate. The heat 
capacity of liquid helium has been discussed on just such 
a basis of elementary excitations and to describe it below 
the critical temperature Landau suggested two classes of 
elementary excitations, phonons or qua nta of longitudi- 
nal compressional waves and rotons [28[. The latter are 
still not completely well understood. The consideration 
of non-static shear phonon contributions to the energy 
spectrum of liquid helium has been largely ignored sim- 
ply because it was not clear whether liquids are actually 
capable of supporting transverse or shear modes. But 
according to Frenkel's proposition a liquid should sup- 
port transverse modes provided the frequency satisfies 
(frequency u > ^, see below). 

Nevertheless the idea that at least longitudinal 
phonons could be excited below 0.5 K has traditionally 
been taken as the explanation for the observation that 
the heat capacity of liquid helium varies as T 3 [2^. In 
most cases investigations have been limited to narrow 
ranges of temperature and pressure these being of imme- 
diate concern in early experiments. Therefore, it is some 
interest to examine the thermodynamic properties of liq- 
uid helium both at elevated pressures, but still presenting 
a liquid phase, and also for wider ranges of temperature. 

Accordingly in this paper we present the heat capac- 
ity of bosonic liquid helium as determined within in the 
framework of a phonon theory of liquids. The physical 
picture of elementary excitations is clarified by means of a 
study of phonon contributions, both longitudinal-like and 
transverse-like, to the energy spectrum of liquid helium 
at elevated pressures. From the analysis of experimental 
data of viscosity and heat capacity [30( and theoretical 
calculations we find that liquid helium remains quantum 
at elevated pressures for a broad temperature range. 



THE PHONON THEORY OF LIQUIDS 

The fact that the solid-like value of heat capacity 
in liquids at the melting point may be summarized by 
Cy = iNks, and the fundamental observation that liq- 
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FIG. 1. Experimental and calculated c„, experimental and 
interpolated rj (insets) for liquid helium. Experimental cv 
and r\ are taken from the NIST database at pressures 20 MPa 
(Fig. [TJa) and 45 MPa (Fig. [TJb). Values of t d used in the 
calculation are 1.41 ps ( 4 He at 20 MPa) and 1.40 ps ( 4 He at 
45 MPa). Values of are 0.0070 GPa and 0.0074 GPa. 
Experimental values of a calculated from the NIST database 
at the corresponding pressures above are 8.7 • 10~ 3 KT 1 and 
8.3- 10 -3 K _1 . Values of a used in the calculation are 3- 10 -3 
K _1 and 2.8T0 -3 K _1 . The uncertainty of both experimental 
heat capacities and viscosities is about 5-10% [3(J. Insets also 
show viscosity fits. 



uids retain the property of fluidity leaves us with an ap- 
parent contradiction. In order to reconcile it, and as 
noted above, J. Frenkel introduced the average time be- 
tween two consecutive atomic jumps thus providing a 
microscopic description of Maxwell's phenomenological 
visco-elastic theory of liquid flow 31]. If r is large com- 
pared with the period of atomic vibrations, a liquid is 
characterized by vibrational states as in a solid (solid 
glass) including shear modes with frequency uj > -. If 
r is small compared with time during which an external 
force acts on a liquid, usually liquids flow. The solid-like 
ability of liquids to sustain high-frequency propagating 
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modes down to wavelengths on the atomic scale, at the 
temperature around and above melting point, was ob- 
served fairly recently 32-36j|. 

As noted, t is a fundamental flow property of a liquid, 
and is directly related to liquid viscosity r/: 77 = G^r 
17|,|3l|, where G m is the instantaneous shear modulus. 
According to time scale the motion of an atom in a liq- 
uid can be viewed as of two types: quasi-harmonic vi- 
brational motion around an equilibrium position as in a 
solid glass, with Debye vibration period of about td =0.1 
ps, and diffusive motion between two neighboring posi- 
tions, where typical diffusion distances exceed vibrational 
distances by about a factor of ten. When r significantly 
exceeds td, the number of diffusing atoms and, therefore, 
the diffusing energy, becomes small and can be ignored. 

The phonon theory of liquids allows us to calculate 
liquid internal energy in general form which can be com- 
pactly presented as 



E = NT 



where D(x) = -% J 


is the thermal expansion coefficient, N is the number of 
phonon states, T is the temperature, and u>d and ujf are 
Debye and Frenkel frequencies correspondingly. 

Eq.([T]) accounts for longitudinal and also for two high- 
frequency shear modes with frequency ui > — . It orig- 
inates at the same level of approximation as Debye's 
phonon theory of solids by using the quadratic density 
of states. The result for a harmonic solid follows from 
Eq.([T]) when ujp = 0, corresponding to infinite relaxation 
time, and thermal expansion coefficient a — 0. This the- 
ory of liquids incorporates the effects of anharmonicity 
and thermal expansion, which is very important not only 
for classical liquids such as Hg and Rb [38J , but also for 
liquid helium as we can see further. 

The phonon theory of liquids has recently been for- 
mulated in a form that predicts the heat capacity of 21 
different liquids, among those: noble, metallic, molecular 
and hydrogen-bonded network liquids. The theory cov- 
ers both the classical and quantum regimes and agrees 
with experiment over a wide range of temperatures and 
pressures [39j]. 



HEAT CAPACITY OF LIQUID HELIUM 

Accordingly we now take a derivative of energy E 
(Eq.dTJ) with respect to temperature T at constant vol- 
ume and compare it to experimental data of heat capacity 
per atom: cy — js^- We have used the National Insti- 
tute of Standards and Technology (NIST) database 
that contains data for many chemical and physical quan- 
tities, including cy for liquid helium. We aimed to check 



our theoretical predictions in a wider range of tempera- 
ture, and therefore selected the data at pressures signif- 
icantly exceeding the critical temperature and pressure 
of liquid helium (T c = 5.1953 K and P c = 0.22746 MPa) 
where it exists in a liquid form in the broader tempera- 
ture range. As a result, the temperature range in which 
we calculate cy is about 40-45 K. Viscosity data was 
taken from the same database (30j | . and fitted in order 
to use in Eq.([T]) to calculate cy. We used the Vogel- 
Fulcher-Tammann (VFT) expression to fit the viscos- 
ity data: r\ = rjg expA/(T — To). To calculate cy from 
Eq.([T]), we have taken viscosity data at the same pres- 
sures as cy and converted it to r using the Maxwell re- 
lationship t = -Q-, where the Frenkel frequency can be 

conveniently expressed as u)p = = 2 "^" ■ 

Eq.([T|) has no fitting parameters, because the parame- 
ters u)£>, a and Goo are fixed by system properties. Val- 
ues of these parameters used in Eq.([T|) are in a good 
agreement with typical experimental values. There is a 
difference between the experimental a and the a used 
in the calculation. At each pressure the experimental a 
was estimated from the formula a = y^p- Experimen- 
tally, VocT only approximately. Consequently, we ap- 
proximated V=V(T) by a linear dependence (an approx- 
imation results in somewhat different a used in Eq.([T])). 
Further, td = used in the calculation (see the cap- 
tion in Fig. [TJ is consistent with the known values for low 
temperature liquid helium under pressure that are typi- 
cally in the 1-2 ps range [4(|. The uncertainty of both 
experimental heat capacities and viscosities is about 5- 
10% 301. 



DISCUSSIONS AND CONCLUSIONS 

As noted earlier, there are two basic analytical ap- 
proaches to calculate liquid energy and heat capacity: 
from the gas phase and from the solid. The approach 
from the classical gas phase has two main contributions 
to liquid energy; kinetic and potential parts and can be 
presented as 



E = K+ gUdV 



(2) 



where K is kinetic energy, g is normalized correlation 
function and U is the interatomic energy. Generally the 
expression in Eq.([2]) is difficult to evaluate for a many- 
body systems and it is not clear how to rigorously in- 
corporate quantum effects at elevated temperatures (say 
tens of Kelvins for He) into Eq. <j2j) . To describe the be- 
haviour of cy on the basis of Bose-Einstein condensation 
(BEC) in liquid helium at low temperatures was origi- 
nally initiated by F. London 4l|. The heat capacity of 
BEC varies as T 3 / 2 at low temperatures and must pass 
through a maximum. This maximum has the character 
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of a cusp which appears at critical temperature T c and 
it has nothing to do with a subsidiary maximum which 
liquid helium possesses at elevated pressures and higher 
temperatures (see Fig. [1]). 

The odd behaviour of the heat capacity of liquid he- 
lium at elevated pressures can now be explained in the 
framework of the phonon theory of liquids. When we 
calculate liquid energy and heat capacity in this theory 
the problem of strong interactions is avoided from outset 
and based on displacive physics associated with phonon 
contributions. We predict that transverse waves exist in 
liquid helium at high pressures. This prediction can be 
verified in a future experimental work. The experimental 
data for heat capacity and viscosity of liquid helium con- 
firms our hypothesis [30j j . Even at elevated pressures liq- 
uid helium persists at temperatures which are 3-5 times 
lower than Debye's temperature (see FigUJ. As the tem- 
perature is raised in liquid helium, more longitudinal and 
transverse- like phonons become progressively excited and 
therefore the heat capacity rapidly grows, which is very 
abnormal for ordinary liquids [39(. Thus liquid helium 
at this P-T region persists as a quantum liquid and also 
as a solid-like liquid with non-static shear rigidity, simi- 
lar to classical liquids. When the Debye and the Frenkel 
temperatures become roughly comparable, cy of liquid 
helium enters the saturation region ('hump'). Further 
increase of temperature then leads to the dissipation of 
transverse-like waves and cy of liquid helium becomes 
shallow (see Fig|T]), implying that at very high tempera- 
tures cy reaches its asymptotic value 3/2. 

The results just presented appear to be fairly accurate 
over the temperature range of experimental importance, 
despite the fact that the formal expression for liquid in- 
ternal energy is quite trivial (see Eq.([T])). The agreement 
is somewhat worse at intermediate part (the slight maxi- 
mum or 'hump' region mentioned earlier) of the cy curve, 
however the maximal difference between the predicted 
and experimental values is actually comparable to the 
experimental uncertainty of cy, namely of 0.1-0.2 J/K 
[30( . But this can also be attributed to an oversimplified 
form of the Debye spectrum of phonon-like states used in 
our analysis. 

Liquid helium in its normal state and at atmospheric 
pressure is actually above the conditions required for 
a Frenkel line (see the diagram 42j) and barely sup- 
ports transverse-like elementary excitations. At signif- 
icant pressures (around 10 GPa) liquid helium quite re- 
sembles the other noble-gas liquids; there are transverse- 
like excitations but at its melting temperature almost all 
phonons are already excited. This state of liquid helium 
is rigid but not quantum. Here we are suggesting an in- 
teresting intermediate pressure range of 5 MPa-500 MPa 
where liquid helium near melting temperature is already 
solid-like but is still significantly quantum. 
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